Photoinduced transition between conventional and topological insulators in 

two-dimensional electronic systems 
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Manipulating the topological properties of insulators, encoded in invariants such as the Chern 
number and its generalizations, is now a major issue for realizing novel charge/spin responses in 
electron systems. We propose that a simple optical means -subjecting to a driving laser field with 
circular polarization, can be fruitfully incorporated to this end. Taking as a prototypical example 
the two-band insulator first considered by Haldane, we show how the electron system can be tuned 
through phases associated with different Chern numbers as the laser intensity is adiabatically swept; 
i.e., a photo-induced analog of the quantum Hall plateau transition. The implications of our findings 
includes the possibility of laser-tuning a conventional insulator into a quantum spin Hall system. 

PACS numbers: 78.67.-n, 73.43.Nq 



The past few decades have seen the topology-related 
quantum aspects of electron systems rise to the forefront 
of condensed matter physics [1]. The unifying concept 
behind these developments is the Berry phase, i.e., an 
extra phase which the wavefunction picks up as the sys- 
tem undergoes an adiabatic excursion. In a semiclassi- 
cal description of Bloch electrons, the relevant excursion 
takes place in momentum space, motivating the intro- 
duction of the Berry connection .4(k) = i(w(k)|Vk|w(k)), 
where u(k) is the periodic part of the Bloch wavefunc- 
tion; a contour integral of .4(k) gives the Berry phase 
associated with that path. Also playing a crucial role 
is the Berry curvature 17 (k) = Vk x w4(k), which adds 
an anomalous- velocity term to the standard equation of 
motion of an electron[2]. The incorporation of these new 
entities has deeply enriched modern solid state physics 
with e.g., a fully quantum mechanical theory of elec- 
tron polarization[3, 4], the notion of adiabatic quantum 
pumping [5], an intrinsic mechanism of the anomalous 
Hall effect [6], and schemes of generating a dissipation- 
less spin current [7]. 

Largely triggered by the discovery of the quantum spin 
Hall effect (QSHE), these activities have in more recent 
years culminated in the conception, classification, and 
an intensive search for new states of matter collectively 
referred to as topological insulators (TI). In its broad- 
est sense, TIs are bulk insulators which nevertheless ex- 
hibit novel linear responses to external fields. These re- 
sponses are dominated by transport through quantized 
edge/surface channels, whose robustness is "protected" 
by a set of topological invariants characteristic to that 
system. A well known example of the latter is the Chern 
number of a 2d system, C = ^ Jg^ (I^Fp'^C^)' where 
the integration is to be performed with respect to the 
first Brillouin zone. Introduced into electron physics by 
Thouless et al in their pioneering work on the integer 
quantum Hall effect (IQHE), this integer- valued quantity 
basically measures the effective "magnetic flux" intrinsic 
to the Bloch electrons, thus giving rise to a quantized 



Hall conductivity a^y = X^n^" !?' ^i^h the summation 
taken over all occupied bands [8]. Hence a Chern insu- 
lator (CI), i.e., a 2d TI with Ccft = E„C„ 7^ 0, real- 
ized in the absence of a magnetic field exhibits a spon- 
taneous quantized Hall effect. States corresponding to 
different values of Ccff are to be regarded as belonging 
to distinct phases, as is the case with each Uxy-plateau 
phase of the IQHE. Haldane was the first to devise a 
minimal two-band model without a magnetic field which 
accommodates CI phases [9]. Despite the artificial na- 
ture of its construction, this model bears immediate rel- 
evance to numerous physical systems (see later discus- 
sions). Furthermore the essence of QSHE is understood 
in terms of a superposition of two copies of Haldane mod- 
els (one for each spin component) with a net time-reversal 
symmetry[10, 11]; recent work clarifies how this picture 
relates to a more generic characterization of QSHE em- 
ploying Z2-valued invariants. 

From the perspective of material function engineering, 
a major shortcoming which plagues the search for TIs is 
that their defining topological invariants, such as Cos, are 
usually material constants whose values cannot be var- 
ied freely within a given sample. This is to be contrasted 
with the situation in IQHE [8], where an adiabatic sweep- 
ing of the magnetic field strength drives the system from 
a normal (Ccff = 0) to a topological (Ccff 7^ 0) insulator, 
and vice versa. In this Letter, we pose and answer in 
the affirmative the question: is there a generic approach 
allowing us to externally tune the topological invariant 
of a TI? Taking up the most basic example of the Hal- 
dane model, we show below that an optical means -the 
application of a monochromatic driving laser field with 
circular polarization, can be invoked to obtain a photo- 
induced normal- Chern insulator transition; i.e., by slowly 
changing the driving laser amplitude (or the laser fre- 
quency), one can transport the electron system into dif- 
ferent regimes of the Ccff vs model parameter phase dia- 
gram. This scenario naturally lends itself to more generic 
topological phase transitions as we later mention, and 
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thereby suggests a route to manipulating the topological 
properties of various insulators [12]. 

Let us then illustrate how the claimed result comes 
about. We begin by noting that our purpose requires us 
to carefully trace the low-energy sector of the electron 
system through its evolution, as the laser field intensity 
is adiabatically swept. This we accomplish by adopting 
a systematic method previously developed by one of the 
authors [13]. Aided by the powerful machinery of the Flo- 
quet theory, the procedure sketched below is applicable to 
the adiabatic dynamics of any tight-binding model, and 
was previously used to give a reliable prediction on how 
the electric polarization of a (normal) insulator renor- 
malizes substantially when subjected to a slowly varying 
laser field[13]. While the polarization is also a prominent 
Berry phase effect, it is worth stressing that the adia- 
batic flow of a global topological quantity such as CcS is 
considerably harder to foresee. Consider then the Hamil- 
tonian for the Haldane model Hh coupled to a circularly 
polarized driving laser field H{t) = Hh + H'{t), where 

HH=tiY.Y. [a\r)b{r + d^) + H.C. 
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The electrons reside on a 2d honeycomb lattice com- 
prising two triangular sublattices which we denote as 
A and B; we associate with each the electron creation 
operators and b'^ . The vector dj connects near- 
est (next-nearest) neighbors [see Fig. Ij. The system is 
at half-filling, and is thus a band insulator. The stag- 
gered potential A and the phase </> are each responsi- 
ble for breaking inversion and time-reversal symmetries, 
which causes the lower-band Chern number to assume 
the set of possible values C = 0,±1[9]. Turning to the 
interaction term H' , we choose the wave vector of the 
laser field to be normal to the lattice plane. The driv- 
ing laser field with a fixed frequency ui can be expressed 
as E{t) = £'r(cosa;<, T sinwi), where r = +1( — 1) cor- 
responds to the left (right) circular polarization. Now 
let the time-dependence of the amplitude Et be weak 
enough to justify the use of the adiabatic approxima- 
tion, which is always possible by preparing a suitable 
laser source. Then the sinusoidally time dependent prob- 
lem conveniently maps onto a time-independent eigen- 
value problem Hf\'^) = e|$) by virtue of the Floquet 
theory [14, 15]; the eigenvalue e is the quasi-energy in the 
language of the latter framework. A straightforward cal- 
culation shows that the Floquet Hamiltonian Hp consists 




FIG. 1: The lattice structure considered in this study. Solid 
(open) dots denote sites belonging to the A (B) sublattice. 
Electrons residing on the A (B) sublattice suffer an on-site 
electron potential of +A (—A); this staggered nature breaks 
the inversion symmetry. The electron hopping energies along 
di and £{ are ti and t2e"^, respectively. 



of 2 X 2 block submatriccs. i.e., Hp = where the 

entries for each block (in the Fourier-transformed bases, 
ak and bk) are: 



(A + mhij)S^n + t2JN{V3X) e'" 
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with Jtv the Bessel function of the N = n — m -th order 
and A = eaET/hw (a: lattice constant). The remaining 



matrix elements, H. 



(mn) 



and H. 



(mn) 



are identical in form 



to Eqs. (3) and (4), respectively, except for the insertion 
of -1 in front of A and (j) in the former, and A and k in the 
latter. This Hamiltonian operates on wavefunctions de- 
fined in the space compositely spanned by the electronic 
degrees of freedom and photons with energy hiu. The di- 
agonal block of the Hamiltonian, is the n-photon 
sector, i.e., the subspace with n photons. 

Now suppose that the condition W <^ hjo Eq is 
met, where W is the electronic band width and Eq = 2A 
the energy gap. In this case the admixture of neigh- 
boring photon sectors is negligible, which enables us 
to concentrate on the zero-photon sector [16], = 

T.k{albl)n{ll 



H = 2i2J^o(V3A)cos0^cos(fc • /j)I 

j 

+ tiJo{X) ^ |cos(fc • di)ax ~ sin(fc • di)<yy^ 



A - 2^2 Jo(\/3A) sin<?!)^ sin(fc • £j) 



(5) 



where 1 is a 2 x 2 unit matrix and Ox^y^z are the Pauli 
matrices. It is straightforward to see that H. is identical 
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in form to the Haldane model, save for the fact that the 
hopping integrals ti and t2 are modulated by factors in- 
volving Bessel functions. Hence, expanding around 
the two Dirac points Ka=± = (a47r/3\/3a, 0), wc finally 
arrive at our effective Hamiltonian 

Ha = -3t2Ja{V3X)cos(f)l + ^tiJo{X) [-ak^a,^ + kyUy] 

(6) 



sin m <7z^ 



where wave numbers k^ y are now measured relative to 
the K± points. Note that the zero-photon sector is in- 
sensitive to the direction of the circular polarization. It 
is an easy task to read off from the results of Ref.9 that 
the Chern number for our effective system described by 
Eq. (6) is given by 



C 



a sgn (A + a3\/3t2^o(\/3A) sin( 



(7) 



This expression for the photo-induced modification of the 
Chern number represents a principal result of this Let- 
ter: in addition to the two variables and A/t2 gov- 
erning the topological property of the original Haldane 
model[9], wc now have a tunable third variable X at our 
disposal which comes from the coupling to the laser field. 
The extended Chern number phase-diagram within our 
three-dimensional parameter space is depicted in Fig. 2. 
Points belonging to the interior of the ravioli sheet-like 
structure fall within one or another of the Chern insula- 
tor phases, with either C = +IorC = — 1. Elsewhere the 
system is a normal insulator with C ~ 0. The array of 
C = +1 and C ~ —1 phases constitute a checkerboard- 
like pattern. By simply sweeping the laser amplitude, 
one can traverse the phase diagram along the A direction 
and cross phase boundaries, thus inducing normal-Chern 
insulator transitions. We have confirmed that the gap 
closes and reopens in the vicinity of the phase transition 
point, as it must, since C is robust to a continuous de- 
formation of the energy bands [12]. We note in passing 
that the circular polarization was essential in the above; 
a replacement with a linearly polarized field would fail to 
yield a massive Dirac-fermion-like energy spectrum as in 
Eq.(6) [17]. 

We now run a consistency check among several con- 
ditions which were implicit above; they will place some 
restrictions on the parameter values for which the tran- 
sition is expected to be observable. We first note that 
|A/t2| < 3a/3, which derives from the observation that 
A = occupies the largest portion of the phase dia- 
gram with C ^ 0. This needs to be consistent with 
the aforementioned condition for the hopping energy, 
W ^ ti hu! Eq- That, however is not straightfor- 
ward, because on top of these \t2Jo{V5^) /tiJo{X)\ < 1/3 
should also be satisfied in order to prevent the bands from 
overlapping [9]. To confront with this problem, wc will re- 
lax the condition W Huj Eq to W < huj < Eq, 
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FIG. 2: (color online). The three-dimensional phase diagram 
of the normal and Chern insulators. (See text for details.) 
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FIG. 3: (color online), (a) Partial section of the phase di- 
agram for A = 0. In the fiUed-area, one can traverse be- 
tween the interior and exterior of the double-sheet structure 
by changing A within the interval < A < A. (b) Section of 
the phase diagram for — ±7r/2. (c) Corresponding Chern 
number as a function of A G [A''', A-^] for <^ = 7r/2 (thick solid 
line) and (f> = —n/2 (thick dashed line). 



which in turn will generally impose a restriction on 
the range of values, [0,A], which the parameter A can 
assume without breaking down our scheme. (Note that 
when W <^ hu! <^ Eq is satisfied, there are no such re- 
strictions on A.) We have evaluated this range by numer- 
ically checking the agreement between the quasi-energy 
bands of the Floquet Hamiltonian terminated at a cer- 
tain order and the eigenenergy bands of When 
t2/ti = 1/3 and A/ti = 1.3, we find that our scheme is 
valid for A ~ 0.3. We therefore confine our discussion to 
this regime in subsequent discussions. The region where 
the normal-Chern insulator transition occurs is displayed 
in Fig. 3(a); in the fiUcd-area, one can traverse between 
the interior and the exterior of the sheet structure by 
changing A within the interval < A < A. 

The normal-to-Chern insulator transition that we 
found can be viewed as a photo-induced analog of the 
inter-plateau transition in IQHE, wherein the laser field 
intensity replaces the role of the external magnetic field. 
(Recall that each IQHE Hall-platcau has a different value 



4 



of Coff.) To strengthen this analogy it is instructive to 
consider sweeping the system between the two points 
A and B in Fig. 3(b), which are the sections of Fig. 2 
at (j) = ±7r/2, and exphcitly display how the Chern 
number change actually occurs. The path connecting 
the two points crosses the phase boundary at a critical 
point A = A*, where a "inter-plateau transition" takes 
place as depicted in Fig. 3 (c) in thick solid (dashed) 
lines. The field strength for the transition around the 
value of A* ^ 0(0.1) is estimated to be of the order of 
^ 0(10^) [V/cm]. This is to be compared, e.g., with the 
corresponding estimate of ^ 0(10^) [V/cm] [18] made 
in relation to dynamic localization, a prototypical exam- 
ple of optical manipulation of electron dynamics. The 
above parallelism loosely conforms with earlier work (not 
straightforwardly applicable to our adiabatic process) 
suggesting that a circularly polarized light should affect 
electron motion in an manner analogous to a magnetic 
field [19]. 

Let us now size up what the implications of our 
strategy are for various physically motivated systems. 

(1) It has been shown[20] that the distribution of the 
Aharanov-Bohm-like flux inherent in the Haldane model 
(described by the phase (p) can be mimicked by tight- 
binding electrons on a Kagome lattice moving in the 
background of a spin chirality ordering. This idea has 
lead to a detailed analysis on the anomalous Hall effect in 
several pyrochlore ferromagnets, such as Nd2Mn207[21]. 
Our method readily carries over to this system; in- 
deed, in the notations of Ref.20 we find that applying 
a circular polarized light to this particular model yields 
C — — sgn ( Jo(A) sin (f)) as the counterpart to Eq.(7), from 
which we expect to encounter the photo-induced transi- 
tion C = — 1— >C = +lasa function of laser intensity. 

(2) It has also been proposed that the Haldane model 
can be simulated by a fermionic cold-atom on an optical 
lattice [22]. A simple estimation reveals that the energy 
scale of the laser fields which constitute the lattice (~ 1 
[eV]) and that of the additional sweeping field 10^^° 
[eV] or larger) are easily separated, which makes the im- 
plementation of the adiabatic scheme feasible. (3) The 
canonical Kane-Mele model[10] for the QSHE, in the 
absence of a Rashba term, reduces to two copies (up- 
spin and down-spin) of Haldane models. (The Luttinger 
model for light /heavy hole bands in a semiconductor also 
reduces to the same model, provided we replace spins by 
a conserved pseudospin quantum number[ll, 23].) The 
sum of the Chern numbers C-j- and for each spin com- 
ponent cancels in accord with time-reversal symmetry, 
while the difference — Cj. needs not. The preceding 
generalizes trivially to this situation, which now leads to 
a photo-induced normal to quantum spin Hall insulator 
(QSHI) transition. When the Rashba coupling is present, 
the Z2-invariant which takes over the role of C can be 
evaluated following Ref.24, using a matrix-valued gener- 
alization of the Berry connection. Here again we find 



that the system can be laser-tuned through a normal- 
QSHI transition. Details on work along this line will be 
reported elsewhere. 

To summarize, we have demonstrated how a photo- 
induced normal-Chern transition occurs in the Haldane 
model. The coupling to a driving laser field with circu- 
lar polarization turns the transfer integrals into tunable 
variables which can result in a change in the Chern num- 
ber. The strategy is reasonably general, and suggests an 
optics-based scheme for manipulating the intrinsic topol- 
ogy of insulators. 
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